A wide band gap metal-semiconductor-metal nanostructure made entirely from 

graphene 
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A blueprint for producing scalable digital graphene electronics has remained elusive. 
Current methods to produce semiconducting-metallic graphene networks all suffer from 
either stringent lithographic demands that prevent reproducibility, process-induced 
disorder in the graphene, or scalability issues. Using angle resolved photoemission, 
we have discovered a unique one-dimensional metallic-semiconducting-metallic junc- 
tion made entirely from graphene, and produced without chemical functionalization or 
finite size patterning. The junction is produced by taking advantage of the inherent, 
atomically ordered, substrate-graphene interaction when it is grown on SiC, in this 
case when graphene is forced to grow over patterned SiC steps. This scalable bottom- 
up approach allows us to produce a semiconducting graphene strip whose width is 
precisely defined within a few graphene lattice constants, a level of precision entirely 
outside modern lithographic limits. The architecture demonstrated in this work is so 
robust that variations in the average electronic band structure of thousands of these 
patterned ribbons have little variation over length scales tens of microns long. The 
semiconducting graphene has a topologically defined few nanometer wide region with 
an energy gap greater than 0.5 eV in an otherwise continuous metallic graphene sheet. 
This work demonstrates how the graphene-substrate interaction can be used as a pow- 
erful tool to scalably modify graphene's electronic structure and opens a new direction 
in graphene electronics research. 
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Patterning a flat graphene sheet to alter its electronic 
structure was envisaged to be the foundation of graphene 
electronics^ The early focus was to open a finite-size 
gap in lithographically patterned nanoribbons, a neces- 
sary step for digital electronics While early trans- 
port measurements supported this possibility,^ it soon 
became apparent that these "transport gaps" originated 
from a series of mismatched-level quantum dots caused 
by the inability of current lithographically to produce suf- 
ficiently narrow, well ordered, and crystallography define 
graphene edges.—"— A working solution to the graphene 
" gap problem" has yet to be formulated, let alone demon- 
strated. We show that in fact such a solution exists, not 
by patterning graphene, but instead by controlling the 
graphene-substrate geometry. 

We have been able to construct a unique, reproducible, 
and scalable semiconducting graphene ribbon with a gap 
larger than 0.5 eV. Using pre-patterned SiC trenches to 
force graphene to bend between a high symmetry (0001) 
face to a low symmetry (112?t,) facet, we produce a nar- 
row curved graphene bend with localized strain. This 
" topologically-deflned" ribbon is a wide-gap graphene 
semiconductor strip a few lattice constants wide that ex- 
tends hundreds of microns long. The strip is connected 
seamlessly to metallic graphene sheets on both of its 
sides. One metallic sheet is n-doped and the other p- 



doped. From this simple morphology, we have not only 
produced a gap suitable for room temperature electron- 
ics, we have demonstrated how narrow Schottky barriers 
can be fabricated in a device-scalable architecture. 

Nano-patterened graphene ribbon arrays 

Epitaxial graphene on SiC is critical for these studies be- 
cause it allows atomic control of where graphene grows 
and what crystallographic orientation the graphene has 
relative to patterned features on the SiC. On the 
SiC(OOOl) surface, graphene always grows rotated 30° 
relative to the SiC < 1010 > direction^ii Using this direc- 
tionality, SiC trenches where etched so that the growing 
graphene had its "armchair" edge parallel to the pat- 
terned steps. This allows the graphene Brillouin zone 
(BZ) orientation to be accurately specified with respect 
to the substrate so that ARPES can measure graphene's 
linear 7r-bands (Dirac cones) on both the (0001) sur- 
face and on facets formed from the sides of the pattern 
trenches. 

During growth, the SiC trench sides facet into a sta- 
ble new surface (see below). Graphene is known to flow 
from the (0001) surface, over the step, and onto these 
facets, ■'^^ We have set the growth time and temperature 
to produce a single graphene layer above the "buffer" 
graphene-like sheet on the Si-faceii^ Both the buffer and 
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FIG. 1. (a) Schematic side view of the side wall ribbons containing two structural graphene sheets (the lower sheet in red is 
referred to as the "buffer" layer) . Grey regions are curved parts of the graphene that are semiconducting. The schematic shows 
the detector orientation to reach both the AT-point of the flat surface and the facet with normal n^. Sample ky and kx{0,ni) 
cuts through the if-points of the (0001) and facet surface show the Dirac cones and constant energy contours, (b) AFM top 
view of side wall grown ribbons showing the long range order, (c) Perspective AFM view of 18 nm deep graphitized trenches 
with a trench pitch of 400 nm. (d) The graphene BZ orientation relative to the armchair edge ribbons. 



monolayer graphene on the Si-face grow over the trench 
edge onto the sidewall facet and then bend back onto 
the (0001) face of the trench floor ji? Figure [Ija) shows 
a schematic of the experimental ribbon geometry. The 
graphene that grows on the facets are referred to as side- 
wall graphene ribbons. An example of part of the post- 
graphitized arrays is shown in Figs. [T] (b) and (c). The 
ribbons are extremely straight over many microns with a 
trench-height RMS roughness less than one 4H-SiC unit 
cell (1 nm). The estimated width of a sidewall ribbon 
is W = d/ s'mAdi, where d is the trench depth and A6i 
is the angle between the surface normals of the sidewall 
facet and the (0001) face. It is this long range order that 
allows detailed ARPES measurements 

Topological changes in graphene's band structure 

We have used ARPES to measure the topological depen- 
dence of graphene's electronic structure as it flows from 
the (0001) surface to the sidewall facets. We do this by 
orienting the SiC steps edges with the 9 rotation of the 
ARPES detector as shown in Fig. [Ija). In this geome- 
try, with the graphene's armchair edge parallel to the SiC 
step edge, the detector measures ky (perpendicular to the 
TK direction) while 9 scans kx{9, fii) along the TK direc- 
tion, perpendicular to the step edges [see Fig.[Ud)]. This 
experimental setup means that we measure graphene's 
if-point band structure for all orientations the graphene 
may take as it curves from the (0001) surface down over 
the trench walls. If a Dirac cone appears in the detector 
at some 9i, we can identify its facet normal hi from the 
relative away from the (0001) iiT-point angle 9o- 
Figure[5l^a) shows a 9 scan at constant energy from a 12 



nm deep trench array. There is no intensity and therefore 
no cones between angles —36° to —33° and —20° to —12°, 
indicating there are no stable facets in this range or that 
the facet size is extremely small. There are, however, 
four regions where cones appear. Figures [Ub)-(e) show 
the ARPES measured cones from these regions. The K- 
point Dirac cone at 9o = —37° is from the graphene on 
the (0001) flat surface [see Fig.[2i;b)]. This is the typical 
band structure of a single graphene layer above a buffer- 
graphene layer. This graphene layer is n-doped by Ed — 
Ep — —0.43 eV, which is the typical doping of a clean, 
ordered Si-face epitaxial graphene filmJ^ 

At small 9 there are two sets of two Dirac cones with 
6'i = ±9° and ±4.5°. These cones correspond to graphene 
that has grown on (2207), and (1103) facets on one trench 
side and their complimentary (2207), and (1103) facets 
that form on the opposing trench side. These large 
area SiC facets developed and are stabilized during the 
graphitization process. The Dirac cones on these facets 
are slightly p-doped as discussed later. The sidewall 
graphene on these two facets are extremely well ordered. 
The Afcy width of these cones ranges from 0.04A for 
36 nm wide ribbons (the same width measured for an in- 
finite graphene sheet on the (0001) surface) to O.OSA 
for 15 nm wide ribbons. Assuming that the Afc width is 
entirely due to angular variations of the facet surface nor- 
mal caused by step edge meandering, the RMS surface 
normal variation, over 1/xm in the y- and x- direction, is 
A(/) ^ 2° and A9 ^ 1.6°, respectively. Keeping in mind 
that these cones are a measured averaged over > 500 par- 
allel ribbons, the facet stability and long range order is 
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FIG. 2. (a) Quasi-Fermi surface from a 12 nm deep trench array after graphitization with the graphene ribbons' armchair edges 
parallel to the step edge. The plot is a constant energy surface (E — Ef = —0.43 eV) in ky and tilt angle 8 (ky is parallel to 
the graphene armchair edge [see Fig. [ijd)]). Dirac cones from each of the angular regions are shown in (b)-(e). (b) Dirac cone 
from the flat (0001) surface. The spectrum is typical of a 1-layer graphene film above a graphene buffer layer. The electron 
doping is 0.43 eV. (c) Typical distorted cone from the intermediate angular range shown in the grey box in (a), (d) Dirac cone 
from the (2207) orientation, (e) Dirac cone from the (1103) orientation. Lines in (b), (c), and (e) are the TB calculated bands 
for a 1.5 nm wide graphene strip between two doped graphene sheets [see supplemental material]. 



exceptionally high. 

The most important feature in Fig. [D^a) is the transi- 
tion between the flat 2-D graphene sheets on the (0001) 
and facet surfaces. This transition is evident in an angu- 
lar range between = — 33°<6'< —20° where a dramatic 
change occurs in graphene's band structure, switching 
from discrete facet Dirac cones to a continuous distri- 
bution (in 9) of broad diffuse cones. An example cone 
is shown in Fig. ^c). It is clear that this cones has a 
gap and a significantly lower slope (velocity oc dE /dky) 
than the Fermi velocity, vp, measured for both the (0001) 
and facet Dirac cones. Figure [SJa) compares an energy 
distribution curve (EDC) (intensity versus energy at a 
constant ky) through the ii'-points for a transition re- 
gion cone and the (2207) facet cone. The spectra density 
for the transition cone decays rapidly below the valance 
band edge, Ey, at E—Ep^—0.5 eV. This is not a doping 
shift since there is no indication of tt bands above -0.5 
eV. Because the position of Ep in the gap is not known, 
the lower limit on the gap size is 0.5 eV. As we'll now 
show, this gapped region is a manifestation of a nearly 
one-dimensional semiconductor. 

Junction Potential Boundary The gapped graphene 
region between the (0001) surface and the sidewall facet 
is unique because it represents a nearly 1-D graphene 



semiconductor seamlessly connected to two 2-D graphene 
sheets on its sides. The nearly 1-D electronic nature of 
the bent graphene can best be demonstrated by consider- 
ing what a hypothetical ARPES experiment would mea- 
sure in a better known 1-D system, a carbon nanotube 
(CNT). If we had enough photon flux to make an ARPES 
measurement on a single CNT with its axis aligned in the 
y-direction, we would only see dispersion along ky. All 
other fc's are not good quantum numbers because of the 
1-D nature of the CNT. In other words, the APRES spec- 
trum would be the same no matter what 9 we rotated the 
CNT about the y-axis, i.e the ky dispersion will be nearly 
independent of 9. 

If we think of the curved graphene region as a par- 
tial arc of a CNT, then the thought experiment is in 
fact identical to the experiment in this work. The only 
difference is that we use 300-500 self-assembled parallel 
coherent ribbons to give sufficient intensity to make the 
measurement possible. Because the step edges are paral- 
lel to the 9 rotation axis, we rotate around the axis of the 
bent graphene by 9. The results for the "edge-CNT" are 
shown in Fig. ^h). Over a 12° rotation the dispersion 
E{ky) is independent of 9, i.e, E is independent of kx{9). 
This simple experiment demonstrates that the gapped 
graphene is confined to a very narrow region of the film. 
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FIG. 3. (a) Comparison of the EDCs from the flat (2207) 
Dirac cone and a cone from the transition region (a lin- 
ear background has been subtracted for clarity). EDCs are 
through ky—Q for a 24 nm wide facet graphene ribbon. The 
transition region cones shows that the graphene valance band 
is shifted 0.5 eV below Ep. (b) EDCs of the cones in the 
transition region. Normalized scans are for different sample 
rotations, 8. All cones in the transition region have the same 
profile, demonstrating the ID nature of the region, (c) Rela- 
tive intensity of the transition cones as a function of 9. Line is 
a flt to the diffraction limit broadening of 30 eV photoelectron 
from a 1.4 nm wide source. 



Unlike the CNT the curved edge is a finite strip. This 
means that the angular distribution of the ARPES is not 
completely independent of 9. The measured intensity is 
modulated by a ^-dependent envelope function of a finite 
size object. We can use this envelope function to estimate 
the coherent dimension of the semiconducting bend. The 
A9 ~12° angular range of the gapped transition spectra 
can be viewed as the diffraction limit broadening of a 30 
eV kinetic energy photoelectrons electron (A — 0.22 nm) 
from a source with dimension d. Figure [Sjc) shows the 
relative intensity of a cone as a function of 9 in the transi- 
tion region. A fit of the measured I{9) to the broadening 
function sinc^(7r(isin0/A) [see Fig. (Hljc)] gives (i=1.4 nm. 
This width is consistent with transmission electron mi- 
croscopy, TEM, measurements that show graphene bend- 
ing over step edges in 5-10 graphene lattice constants. 
While local probes like STM will be able to measure this 
width more precisely, it is clear that the semiconducting 
graphene strip is very narrow. 



Of course, the bent graphene band structure cannot be 
understood in the context of a CNT since the wave func- 
tion boundary conditions are different in the two system. 
Instead of cyclical boundary conditions in a complete 
CNT, the bent graphene can be considered as a finite 
arc that is continuously bounded by two flat graphene 
sheets on both sides with different doping. While a self 
consistent theory will need to developed, there are several 
important effects that will influence the band structure 
of this system; strain, finite size effects, local coulomb 
potential, and substrate bonding. 

The strain and electronic energies are comparable in 
this system. From Norimatsu and Kusunoki's work, the 
known radius of curvature of a graphene layer over a 
step is ~ Inm The strain associated with this amount 
of bending can be very large. Calculated band-gaps for 
similarly strained graphene can range from 0.2-0.5 eV 
depending on the model.— liii^i Comparing this strain 
energy with a finite size band-gap expected for a 1.4 
nm graphene ribbon. Eg ~ leV-nm/1.4nm = 0.8 eV^^ 
shows that both effects lead to band gaps similar to 
the experimental value. Because the energy scales are 
similar, Peierl's like distortion are potentially important. 
This is an important point because, unlike a CNT or a 
finite width ribbon, the 1-D strip's boundary conditions 
are not discrete. This allows the graphene at the bend 
to flex in order to lower its total energy. 

Besides the effect of strain and boundary conditions, 
there is an onsite potential at the bend due to the dif- 
ferent doping on either side of the ribbon. The charge 
distribution and the Schottky barriers that can form at 
the metal-semiconducting will require a self consistent 
solution to the Poisson equation. This charge distri- 
bution affects the electron boundary condition and will 
be important in determining resonant states that influ- 
ence transport through the metal-semiconductor-metal 
region. 

Because the Schottky barrier is so narrow, tunneling 
through the junction should be an important part of 
graphene transport over steps. Evidence for this, in the 
form of a decrease in graphene's conductivity over ran- 
domly oriented 1 nm SiC steps, has been reportedi^^ Al- 
though graphene's strain going over a 1 nm step should 
be significantly smaller than the taller steps studied in 
this work, the transport findings in Ref. [19| are consis- 
tent with a semiconducting strip at the step. 

While it is obvious that a good theoretical understand- 
ing of the electronic structure remain to be formulated, 
we have used a simple tight binding calculation that re- 
produces many of the salient features of the ARPES 
data [see supplement]. The model consists of a 2000 
atom wide flat graphene sheet. A 1.5 nm strip in the 
middle of the 2D fllm has a different hopping integral 
than the rest of the sheet to mimic the effect of strain 
at the bendingj^S^iS^ In addition, a position dependent 
on-site potential, V{x), is used. This is to account for 
the transition from the n-doped (0001) graphene and the 
p-doped sidewall graphene [see below and in the supple- 
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FIG. 4. (a) A magnified view of the graphene Dirac cone 
from the stable (1103) facet Fig. EJe) near Ef- The figure 
shows EDC offset in intensity to highhght the region of the 
cone near Ep. The facet is p-doped by 70 meV±20 meV. (b) 
EDCs through fcy = of the (1103) facet cones. Data is for 
three different sidewall ribbon widths; 15, 24, and 36 nm. 



mental material]. The band structure calculated using 
this model is shown in Figs.t^Ib), (c), and (e). Note that 
the model produces a significant number of state within 
the cones from the narrow semiconducting region. This 
correlates well with the experimental diffuse intensity ob- 
served within the cone. 

Facet band structure 

The morphology shown in Fig. [IJa) consists of a wide 
band gap, nearly 1-D, semiconducting wire bounded by 
two flat graphene sheets. One sheet on the (0001) surface 
and the other on the sidewall facet. Because boundary 
conditions will be important in understanding the details 
of the ribbon band structur, we must understand the dif- 
ferences between the graphene grown on the (0001) sur- 
face and graphene grown on the sidewall facets. The two 
sheets each are in different charge states. Graphene on 
the (0001) surface is metallic and n-doped. The sidewall 
facet graphene has a small p-doping. Figure SJa) shows 
that the (1103) facet has a ^70 meV p-doping. It is not 
clear from these studies if the narrow sidewall ribbons, 
have a gap, particularly the 15 nm wide ribbons. Figure 
m^b) shows EDCs through the graphene Dirac cone on 
the (1103) facet for different ribbon widths. All ribbons 
show the same shape near Ep. The expected gap for a 



w = 15 nm ribbon is Eg ~ leV-nm/15nm = 66 meVJ^ 
Since the doping is comparable to the expected gap, it 
is difficult to say with certainty whether or not the facet 
graphene is metallic or semiconducting. We can say that 
if there is a gap, it is less than 70 meV. 

Why there is a doping difference between graphene on 
the (0001) and graphene on either the (2207) or (1103) 
SiC surfaces is worth attention because it suggests other 
ways to spatially control graphene's properties. The 
surface reconstructions are different on the two facets. 
LEED shows that while the graphene diffraction spots are 
visible on the stable facet walls, the (6-\/3x6-\/3)R30° re- 
construction spots of the (0001) face disappear [see sup- 
plemental material] . The surface reconstruction must be 
different on the two facets because the bulk terminated 
sidewall facets have a rectangular symmetry with two 
dangling bonds per surface atom while the (0001) sur- 
face has hexagonal symmetry with one dangling bond per 
surface atom. The different surface symmetries implies 
that polarizability, bond re-hybridization, chemical reac- 
tivity, etc., can all be different and influence the charge 
transfer. While thorough studies of these facets remain 
to be done, we can show that the work function for the 
sidewall facet is very different from that on the (0001) 
surface. 

A closeup electric force microscope (EFM) image of 
the ribbons in Fig. [5l^a) shows the work function con- 
trast change on the sidewalls that occurs after graphene 
has grown on the sidewalls. From the photoelectron sec- 
ondary electron cut off, we have measured a ~ 1.5 eV 
decrease in the work function of the facet relative to the 
polar (0001) face. This implies that the buffer graphene 
layer on the (0001) surface reacts very differently than the 
graphene layer above the SiC sidewall facet. This is con- 
firmed in Figs.[Sl[b) and (c) that show a surface-derived^ 
set of bands near the iiT-point that both shift to higher 
binding energy and change shape on the facet. We note 
that the shift down is nearly equal to the work function 
shift even though the Dirac cones remain at the same 
energy. In some sense the reconstruction of the facet, 
with a presumed interaction change with the graphene 
buffer layer, behaves similar to the hydrogen intercalated 
SiC(0001)-graphene system)^ where the hydrogen affects 
the buffer-SiC interface in a way that effectively neutral- 
izes the charge transfer to the graphene film. 

One last point must be mentioned. The bilayer 
graphene pair on the large facets should have an un- 
usual stacking. Assuming that the two sheets bend over 
29° (the angular distance to the first fiat facet), the top 
graphene layer must slip relative to the second layer by 
~ 30% of the graphene unit cell. The strain energy to 
keep the two layers in an AB commensurate stacking with 
this amount of slip is too large and cross-sectional TEM 
shows no evidence for buckling or any other strain reliev- 
ing mechanism that would be required to bring the two 
sheets into commensuration as they go over the step edge. 
It has recently been shown that sliding graphene layers 
can significantly distort the K-po'mt band structurci^ 



6 




Graphene on SiC(OOOl) 




Graphene on SiC(2207) 




-0.4 -0,2 0.2 0.4 



ky (A-') 

FIG. 5. (a) An EFM image of the side wall graphene with 36 
nm wide facets. Light lines are areas were the work function 
changes when graphene has grown on the sidewalls. (b) and 
(c) are valence band spectrum near the Tf-point for the fiat 
(0001) and (2207) facet surfaces, respectively. The dashed 
line shows the surface derived band change. 



This may explain while the experimental two layer Dirac 
cones are similar to those from a single graphene sheet 
and do not show the AB split bands associated with 
Bernal stacking. A closer inspection of the cones from 
the large stable facets does shows a large spectral inten- 
sity ~ 0.5 eV below Ep [see Fig. EJa)] that would be 
characteristic of a disorder induced energy-broadening of 
the split band from an AB stacked pair. However, no sig- 
nificant band broadening in ky, which must be associated 
with the AB band, is observed. 

Conclusion 

We have shown that a nearly 1-D, large band-gap, semi- 
conducting graphene ribbon, seamlessly connected to 
metallic graphene sheets, can be produced from a scal- 
able bottom-up approach by forcing topological changes 



in the SiC that in return affect graphene's electronic 
properties. The ribbons produced this way have a band- 
gap energy > 0.5 eV and are both narrow (^ 1.4 nm) 
and continuous over macroscopic lengths. In a sense, the 
macroscopic arrays of graphene metallic-semiconducting 
junctions presented in this work are the ideal architec- 
ture envisioned for carbon nanotube electronics; perfectly 
aligned ribbons with identical chirality patterned over 
macroscopic distances with atomic precision. The rough 
edge problem associated with lithographically patterning 
graphene ribbons is solved using a the bottom-up assem- 
bly as the graphene grows The "edges" in the struc- 
tures presented here are atomically perfect because they 
are simply a curved distortion in a continuous graphene 
sheet that is defined by potential gradients (elastic or 
electronic) rather than ill defined edge terminations. 

Graphene, and therefore graphene electronics, is too 
often viewed as an isolated entity. In this work we do not 
merely show the importance of the graphene-substrate in- 
teraction, we demonstrate how the support substrate can 
be used to purposely alter graphene's electronic prop- 
erties. We give an important example of how a well- 
ordered, commensurate substrate interaction can be used 
to force topological changes in the graphene that lead 
to metal-semiconducting transitions and alter its doping. 
The large number of ways the substrate can be altered 
means that a new field of scalable geometry dependent 
device research is just beginning. For instance, we have 
only studied armchair edge 1-D wires. Like the chiral- 
ity in CNTs the device properties could be very differ- 
ent for zigzag edge bends. Another example is doping 
graphene. In this work the facet doping is controlled 
by intrinsic SiC surface effects. Because the facet walls 
have a different chemical reactivity, adsorbates can be 
designed with self-directed bonding to sites only present 
on the facet. This means that graphene dopants can be 
patterned with atomic control simply by cutting trenches 
in different crystallographic edges. Exploring the trans- 
port properties of the junctions described in this work is 
just beginning. New architectures based on the princi- 
ples outlined here are sure to significantly alter graphene 
electronics research. 
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Methods 

The substrates used in these studies were n-doped n = 
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2 X 10"'^^cm~^ 4H-SiC. The graphene ribbons arrays were 
made by a selective growth methodJi Vertical trenches 
are first produced in the SiC(OOOl) (Si- face) surface 
and then graphitized in a controlled silicon sublimation 
furnace. The arrays were prepared by first producing 
a negative mask of hydrogen silsesquioxane (HSQ) by 
e-beam lithography. The SiC(OOOl) substrate was ex- 
posed to reactive ion etching with a SFg— O2— Ar plasma. 
This produced 7.5 lun to 18 nm deep trenches depending 
on etching time. The samples were then graphitized by 
heating to 1560 °C in a carbon RF furnace. The high 
density parallel trench arrays were made with a pitch be- 
tween 100-400 nm over a 1 x Smm^ area (2,500-10,000 
1mm long trenches). This allows 300 to 500 ribbons to 



be easily localized in the ~ 40 — 50/im ARPES beam. 

The samples were transported in air before introduc- 
tion into the UHV analysis chamber. Prior to ARPES 
measurements the graphene films were thermally an- 
nealed at 800 °C in UHV. ARPES measurements were 
made at the Cassiopee beamline at the SOLEIL syn- 
chrotron in Gifs/Yvette. The high resolution Cassiopee 
beamline is equipped with a modified Petersen PG- 
monochromator with a resolution E/ IS.E ~ 70000 at 
100 eV and 25000 for lower energies. The detector is 
a ±15° acceptance Scienta R4000 detector with resolu- 
tion A£: < 1 meV and Afc ~ O.OlA"^ at fiuj = 36 eV. 
All measurements were carried out at lOOK. The total 
measured instrument resolution is (Ai?<12 meV). 
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LEED from patterned arrays 

Because the ribbon arrays are highly ordered, it is possible to get meaningful information about the facet structure 
from LEED. Figure [T] shows a LEED pattern from parallel ribbon arrays with 20 nm trench depth. The main features 
of the pattern are the graphene rods from the (0001) surface with the typical 6V3 patteruF^ The graphene is rotated 
30° with respect to the SiC < 0110 > direction. The trench edges are oriented perpendicular to the SiC < 0110 > 
direction so the graphene grows over the trench edges with the graphene < 110 >g direction perpendicular to the 
steps. At the electron energy used in the figure, the SiC spots are reduced in intensity both from a structure factor 
effect and the electron attenuation through the graphene. 

Note the streaking of the (00), (01)g, and (10)g spots. The streaking is typical of a faceted surface.^ It is important 
to realize that only the graphene spots are streaked. This conclusively shows that the graphene is on the sidewalls. 
The GVS spots show no streaking as expected because this pattern is only associated with the (0001) surface graphene. 

Because the position of the specular spot (momentum parallel to the surface is zero) is independent of energy, the 
specular spot from the facet surface can be identified by changing the electron energy and finding the spot whose 
position is independent of energy. From such a series of LEED images, we can identify the specular spot from graphene 
on a SiC facet at an angle of ^29° relative to the (00) rod, i.e, the (2207) facet. This marked in the insert in Fig. [1] 
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FIG. 1. (a) LEED pattern from a graphitized trench array with side wall graphene. The trench depth is 20 nm and the trench 
pitch is 100 nm. The incident energy is 65 eV. The the spots index with subscript "G" are graphene spots oriented 30° relative 
to the SiC(OOOl) surface. The 6\/3 and the (Ol)sic spots from the (0001) surface are also marked. Intensity streaks run through 
the graphene and (00) rods. Inset shows the specular rod from the (2207) facet. 



Tight binding calculations 

An estimate of the band structure of the 1-D bent graphene, continuously connected between two metallic graphene 
sheets with different doping, was calculated using a tight binding method in the nearest-neighbor approximation. 
We model the system using a 1.5 nm flat graphene nano- ribbon, CNR, with armchair edges to represent the bent 
graphene. This strip is connected to a 17.25 nm wide p-doped strip representing the side wall graphene and an n-doped 
104.15 nm wide strip representing graphene on the (0001) surface [See Fig.[5Ja)]. The total number of carbon atoms 
across the three regions is 2000. The model is periodic in the y-direction, perpendicular to the armchair edge. 

For the purposes of demonstrating the confining effect of strain in the bend, we assuming a position dependent 
hoping parameter, t{x)A'— We allow t{x) to be different in the narrow ribbon compared to the large graphene sheets 
on either side. To include the different doping of the three graphene sheets, we allow a spatially varying on-site 
potential Vi{x). The hamiltonian describing the system is then; 



H = -t{x) a+bj + V.,{x)h, + h.c, 



(1) 
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where (a^) and bj are the annihilation (creation) operators for the sublattice A and B, respectively. The 

operator fii is the local charge density operator— and is given by fii — afdi {hi = bfbi) for the A and B sublattice. 
FigureslDJb) and (c) show the spatial variation of t{x) and Vi{x), respectively. Because calculating Vi{x) in the ribbon 
is beyond the scope of this calculation, we describe the potential in the ribbon with a 50 meV deep well, smoothly 
connected to the two differently doped border graphene sheets as shown in Fig. [2l While the shape of this well 
influences the resonant states in the system, it does not change the important differences of the band structure in the 
three regions. 
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FIG. 2. (a) A model of the metallic-semiconducting-metallic graphene ribbon. The p-doped part spans a; = — 17.25 nm, 
curved part x = 17.25 — 18.75 nm and n-doped part x = 18.75 — 122.9 nm. The variation of (b) the hopping parameter t{x) 
and (c) the on-site potential Vi{x) across the GNR. We use 1000-GNR (TV — 1000). The length is in the units of graphene C-C 
bond a — 1.42 A. The hopping parameters in different parts of the ribbon were determined from the measured slope of E{k) 
dispersion. The on-site potentials describe the n- {Vi^„ — —0.43 eV) and p- {Vi^p = -1-0.07 eV) doping, smooth potential variation 
in the curved part of the ribbon and an additional on-site potential [Vc — 50 meV) induced by the curvature R = 3 nwm. 



We solve the eigenvalue problem Hijj^^k^ — E^^k^{ky)^i,,k^ in the armchair nano-ribbon geometry using 2N = 2000 
carbon atoms. This leads to 2000 bands, E{ky), that are indexed by the quantum numbers v and k^- The quantum 
number indexes the sub-bands in the conduction and valence bands, is the wave number of the standing waves in 
the narrow ribbons (ID limit) and becomes the Bloch wave number of the propagating states in the limit of infinitely 
wide ribbons (2D limit where GNR becomes a graphene sheet). 

In order to compare the results of this model to the ARPES data, we have to take into account the selectivity of the 
ARPES with respect to kx [see Fig. 1 in the main text]. To do this, we must project the calculated states ipu,k^ onto 
the states \kx^exp 

>, where kx,exp is the the component of k^ parallel to the local surface normal. kx,exp is selected 
by the angle 9 of the detector with respect to the (0001) facet. This is important for properly representing the band 
structure at the JsT-point of the two flat sheets. This projection is not necessary on the curved part of the ribbon 
because the 1-D character of the ribbon. In that case, all the states V'l'.fcx foi' fhe narrow ribbon are mapped onto the 
ID Brillouin zone given by ky = {—■^, so that the measured band structure is independent of 6. 
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